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Two-temperature accretion flows in magnetic cataclysmic 
variables: Structures of post-shock emission regions and 
X-ray spectroscopy 
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ABSTRACT 

We use a two-temperature hydro-dynamical formulation to determine the tempera- 
ture and density structures of the post-shock accretion flows in magnetic cataclysmic 
variables (mCVs) and calculate the corresponding X-ray spectra. The effects of two- 
temperature flows are significant for systems with a massive white dwarf and a strong 
white-dwarf magnetic field. Our calculations show that two-temperature flows predict 
harder keV spectra than one-temperature flows for the same white-dwarf mass and 
magnetic field. This result is insensitive to whether the electrons and ions have equal 
temperature at the shock but depends on the electron-ion exchange rate, relative to 
the rate of radiative loss along the flow. White-dwarf masses obtained by fitting the 
X-ray spectra of mCVs using hydrodynamic models including the two-temperature 
effects will be lower than those obtained using single-temperature models. The bias is 
more severe for systems with a massive white dwarf. 
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1 INTRODUCTION 



Ishida et alJll99ll: IWu et alJll995t iFuiimoto fc Ishid3ll997t 
Ezuka fc Ishidalll999T) " 



Magnetic cataclysmic variables (mCVs) are close binaries 
containing a magnetic white dwarf accreting materi al from a 
Roche-lobe filling low -mass companion star (see e.g. IWarnerl 
ll994ICroDDerlll99(t . The material flow is supersonic when 
it leaves the inner Lagrangian point of the binary. It be- 
comes subsonic near the white-dwarf surface, and a shock 
is formed, heating up the accreting material to tempera- 
tures T « 3GA/ w m H /8fci? w ~ 10-50 keV (where G is 
the gravitational constant, k the Boltzmann constant, m H 
the hydrogen- atomic mass, M w the white-dwarf mass, and 
-R w the white-dwarf radius). The accreting material in the 
pre-shock flow is thereby photoionised. The post-shock flow 
is cooled by the emission of bremsstrahlung X-rays and cy- 
clotron optical/IR radiation, as the material settles onto the 
white-dwarf atmosphere. The X-ray emission from a mCV 
depends on the temperature and density structures of the 
post-shock region, which in turn depends on the proper- 
ties, mainly the temperature, of the accretion shock. As the 
shock temperature is determined by the mass and radius 
of the accreting white dwarf, we ca n infer the mass of the 
white dwarf from the X-ray spectra iRothschild et aflll98lt 



The X-ray spectra of mCVs, in particular the subclass 
intermediate polars (IPs), are well fitted by model spec - 
tra generated by the Aizu-type shock models lA^zu|^73), 
such as those in IChevalie r fc Imamural i l982T) ; IWul il994l) ; 
IWu et al.1 (11994) and lCropper et alJ <1999l) . In these mod- 
els, the electrons and ions have the same temperature locally 
(here termed a one-temperature model), but the tempera- 
ture and density change along the flow in the post-shock 
region. It has been noticed that the white-dwarf masses ob- 
tained by the X-ray spectral fits using these models tend to 
be systematically larger than those deri ved from some other 
methods, e.g. optical spectroscopy fsee lRamsav et ai]ll998l : 
iRamsavlEoOoT) . 

The discrepancies could be due to inaccuracies in these 
other determinations. Alternatively, they may arise when 
the assumed absorption column density in fitting the X- 
ray data is uncertain. They could also be due to the fact 
that some relevant processes have not been considered in 
deriving the temperature and density structures of the post- 
shock emission region. In the one-temperature model, the 
electrons (the radiating particles) and the ions (the major 
energy-momentum carriers) are strongly thermally coupled 
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and share the same temperatures. The coupling is main- 
tained by electron-ion collisions. However, if the radiative 
cooling timescale of the accretion flow is shorter than the 
electron-ion collision timescale, the electrons will lose their 
energy rapidly while their energy gains via collisions are un- 
able to keep up radiative loss. The two charged species then 
have unequal temperatures, an d the accretion flow becomes 
a two-temperature flow. (See Iwiill2000l: lBeuermannll2004 
for recent reviews of accreting shocks in mCVs.) 

Two-temperature flows are expected to occur in mCVs 
containing a wh ite dwarf with a very st rong mag netic field 
(B > 30 MOilLamb fc Masters! <1979ft fsee also llmamural 
I l98 illmamura et alJl987ill996tlWoelk fc BeuermanrJl996t 
Isaxton fe WJl999l . l200ll:lFischer fc Beuermanr]|200lt) . Cv- 
clotron cooling is efficient in these systems, and two- 
temperature effects are most significant in the down-stream 
region just beneath the accretion shock. The electron tem- 
perature and density structures of a two-temperature ac- 
cretion flow can be very different from those of a one- 
temperature flow. Moreover, the thicknesses of the high- 
density X-ray emitting region in the two flow models are 
expected to differ as well. Thus, the mass estimates ob- 
tained from a one-temperature post-shock model and a two- 
temperature post-shock model could be different. 

Here we investigate the X-ray spectral properties of 
two-temperature accretion flows in mCVs. We consider a 
semi-analytical approch, in which the hydro dynamic models 
are co nstru cted following the pre scriptions of llmamura et ail 
J1987D and lSaxton fc Wul feOOll) . The thermal coupling be- 
tween the ions and electrons is parametrised by the Coulomb 
collision rate. The radiative loss is due to the emission of 
bremsstrahlung X-rays, which is optically thin, and opti- 
cal/IR cyclotron radiation, which could have substantial 
optical depths. The total radiative loss is approximated 
by a composite cooling func tion a s in t h e previous studie s 
of one-temper ature flows by [w] Jl994l) : IWu et alJ Jl994) : 
ISaxtonl ill 9991) . We calculate the temperature and density 
structures of the two-temperature post-shock emission re- 
gions and generate the X-ray line and continuum spectra 
(S by con volving the MEKAL optically thin therma l plas- 
mas model jMewe et alll985llKaastra fc Meweil993l) in the 
XSPEC package. The results of the two-temperature calcula- 
tions are compare d with the results of one-temperature cal- 
culations (see e.g. | Cropper et allll998t iTennant et aljfl998t 
ICropper et alJll999f) . 



2 POST-SHOCK ACCRETION FLOW: A 
TWO-TEMPERATURE FORMULATION 

Our hydrodynamic formulation assumes that the gas in the 
post-shock region is completely ionised. The flow is along 
the magnetic field lines. We omit the gravitational force and 
the curvature of the field lines. These effects could be im- 
portant when the thickness of post-shock region is signifi- 
cant in comparison with the white-dwa rf radius. (See e.g. 
ICropper et aDll999l ICanalle et aT1l2005h Thus, the flow is 
perpendicular to the white-dwarf surface and is practically 
one-dimensional. Futhermore, we consider only the station- 
ary situation. Hence, the hydrodynamic equations governing 



the flow are 

v dp 
dx 

dP c dPj 
dx dx 

dPc 

dx 

dPi 



+ P 



pv 



dv 
dx 
dv 
dx 
vP c dp 
p dx 
vPi dp 



(7 



IT - 7— = -(7 



dx 



p dx 



) (r ci 
1) r c 



•A) 



(i) 

(2) 
(3) 
(4) 



where v is the flow velocity, p the density, P e the electron 
partial pressure, Pi the ion partial pressure, F e i the rate 
of the electron-ion energy exchange, A the electron cooling 
function, and 7 the adiabatic index. The total gas pressure 
P is the sum of the electron and ion partial pressures. 

We assume an ideal gas law for both the electron and 
ion gases, i.e., 7 = 5/3 and P{ c ,i} = ^{ c ,i} kT{ ei q, where n e 
is the electron number density, rii the ion number density, 
T e the electron temperature, and Ti the ion temperature. 
The rate of energy exchange due to electron-ion collision is 
r o i m 3nifc(Ti — T c )/2t i, where t e i is the equi-partition time, 
given by 



tei — 



3m c miC 



(0c + 0i) 



3/2 



(5) 



8 (2-Ky/ 2 Zfn c e 4 In C 

JSpitzerl 11963) . where m e is the electron mass, m Y the ion 
mass, c the speed of light, e the electron charge, 9 C — 
kT e /m c c 2 , 8i — kTi/rriiC 2 , and In C the Coulomb logarithm. 
We have 



(6) 
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The composite cooling function A consists of a 
bremsstrahlung cooling t erm Ab r and a cyclotron cooling 
term A cy dWu et alJll994l) . i.e., 



A = Abr + A c 



A b 



2k 



n c , s 



(7) 



where e = thi/t cy is the ratio of the bremsstrahlung-cooling 
timescale to the cyclotron-cooling timescale. (Here and else- 
where, the subscript s denotes the value evaluated at the 
shock surface.) The derivation of the expression of e B is given 
in Appendix IA3I The explicit form of the bremsstrahlung 
cooling function is 



A br = 16 



2tt\ 



3/2 



r?2 a l/2 

2 Z f n c m 0J g B 



(8) 



m e c 2 h 

jRvbicki fc Lightmanlll979l) . where h is the Planck constant 
and g B the Gaunt factor. 

We note that our calculations depend on the functional 
form of A cy in addition to the parameter e s . In J7J we have 
assumed a power-law type function to approximate the cy- 
clotron radiative loss term. A key ingredient in constructing 
the cooling function is to estimate the fre quency a;* at which 
the local cyc lotron spectrum peaks (see lWada et aljll98(J : 
ISaxtonll999l ). How well the assumed power- law function ap- 
proximates the cyclotron radiative loss depends on the accu- 
racy of determining ui* in a given geometrical a nd hydrody- 
namic setti ng. Her e we fo llow the approach of IWada et alJ 
Jl980h and lSaxton! ill 9991) . but the technique would be im- 
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proved if we could construct a local cyclotron cooling func- 
tion more self-consistently, say by using an iterative scheme 
that calculates the cyclotron emission spectrum of each stra- 
tum and uses it to readjust the local parameters of the cy- 
clotron cooling function. 

We ignore electron conduction, Compton scattering and 
nuclear burning in the energy transport. These processes are 
unimportant in the accretion flows of mCVs, unless the sit- 
uation is very extreme, e.g. the white dwarf is unusually 
massive (~ 1.2— 1.4 Mq) and the accretion rate is very high 
(M > 0.1M P , where M E is the Eddington accretion rate) 
ilmamura et al.ll987l) . We do not include line cooling in our 
calculation of A. However, line cooling may not be negligible 
at the very bottom of the post-shock region where the tem- 
perature is low. For systems with low white-dwarf masses, 
the shock temperature and the post-shock gas temperature 
are low enough that the Fe L lines can actua lly contribut e 
a significant fraction of the total cooling (see Mukai 2003). 
A fully consistent treatment of line cooling in the hydrody- 
namic calculation is non-trivial, and we will leave this for 
future studies. 

To simplify the hydrodynamic equations, we consider 
the dimensionless variables £ = x/x s , r = —v/va, £ = 
p/pe,, 7Ti = Pi/p^Vff, and 7T C = Pe//O a «g, where v e = 
(2GM„ / Rw) 1 / 2 is the free-fall velocity at the white-dwarf 
surface, p a = m/va the density of the pre-shock flow, and m 
the specific accretion rate. Substituting these variables into 
equations (1), (2), (3) and (4) yields 





= o, 


(9) 


d-K c d-Ki dr 

~di + ~di + T ^di 


= o, 


(10) 


dn c T7r c d( 


= A — f ei , 


(11) 


dlTi TTT[ d( 


= f ei , 


(12) 



where we define expressions for the non-dimensional energy 
exchange and cooling functions, r c i = (7 — l)[x a /peVg) r e i 
and A = (7— l)(a; 8 /pa«ff) A. 

The boundary values for electron and ion pressures (7r e , a 
and 7Ti :S ) are determined by the efficiency of the electron-ion 
coupling through the shock transition region. Their ratio, 
<t s = 7r CiS /7ri iS , can take values between m e /m,i (the ratio of 
the electron mass to the ion mass) and pi/p e (the ratio of 
the molecular weight of the ions to that of th e electrons), de- 
pend ing on the assumed coupling processes dlmamura et alJ 
1996). The physics of how the electrons couple with the 
ions at the shock is not well understood. The value of <r s 
is really a dependent property of the pre-shock flow, but its 
derivation would require solving a comprehensive model of 
pre- and post-shock regions, with explicit radiative trans- 
fer and hydrodynamics: a complex task beyond the scope 
of the present work. We therefore treat cr s as a free param- 
eter and consider a few sensible values in our calculations 
to see what difference it will make to the density and tem- 
perature structures of the flows and the associated X-ray 
spectral properties. 

For the other boundary conditions, we assume a strong- 
shock condition: at the shock surface (£ = 1), r s = 1/4, 
Q = 4 and for the (unitless) total gas pressure (ttq = 7r c + 7Ti) 



we have 7ro, s = 3/4. At the bottom of the flow, we consider a 
"stationary wall" condition: r = at the white-dwarf surface 
(£ = 0). We omit the transf er of energy from ben eath the 
white-dwarf atmosphere fsee lWu fe Cropper 2001) 



3 TEMPERATURE AND DENSITY 

STRUCTURES OF THE POST-SHOCK 
REGION 

We integrate the mass continuity and the momentum equa- 
tions, yielding 



no 



= 1, 
= 1 



(13) 
(14) 



Substituting these into the energy equations and eliminating 
7Ti, we obtain two differential equations: 



dr 

d-JTc 

where 



A 



7(1 — r) — r 



A — , 7T7r ° A -f 
7(1 — r) — t 



A = (7 - l)x s pv e 2 A J [1 + e s /(r, tt 

V T° 



(7 - l)x e pa.V B X- 



1/2 3/2 



(15) 



(16) 



(17) 



(18) 



the constants A and X depend on the composition of the 
plasma (see Appendices IA1IA21 . 



AOC + /3 



(19) 



describes the efficiency of the secondary cooling process (e.g. 
cyclotron cooling) relative to thermal bremsstrahlung cool- 
ing. The constant \~ (Z + l)/Z depends on the abundance- 
weighted mean charge of the ions. We define a parameter 
ipd — X/Avfi to express the efficiency of ion-electron energy 
exchange compared to the radiative cooling. 

The velocity, density, temperature and pressure profiles 
of the ions and electrons in the post-shock flow can be ob- 
tained by integrating equations 1151 and 1161 . 

Assuming a neutral balance of electron and ion charges 
within the plasma, the one-temperature condition, T c =71, 
implies a particular ratio of the partial pressures, 7r c = Ztti. 
By setting 7 = 5/3, a s = Z (see equation IA7I in the ap- 
pendices) and 7r e = (1 — t )/x m equation I15L we recover 
equation (2) in|w3 dl994Tl (with a = 2.0 and /3 = 3.85) for 
the one-temperature flows with a power-law approximation 
to cyclotron cooling. If we assume the same expression for V, 
CT a and 7r c in equation 116H . we will obtain an equation differ- 
ing from equation 1151 by a factor of 1/x in the terms at the 
left side. This is because the assumptions of charge neutral- 
ity and equi-partition between the electron and ion energy 
(i.e. ti c = Z-k\) in the entire shock-heated region requires 
the exchange term to be determined by the electron cool- 
ing rate. An additional assumption of the energy-exchange 
rate depending on the difference between the ions' and the 
electrons' temperatures will inevitably lead to inconsistency. 
If we replace the energy-exchange term in equation H6H by 
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Figure 1. Stationary structures of accretion shocks with (o"s,V'ci) = (Z, 10) and solar metallicity (Z = Zq). The parameter that sets 
the relative efficiency of cycloton cooling to bremsstrahlung cooling is e B = 0, 1, 100, from top to bottom. The first column shows the flow 
velocity t. The second column shows the electron pressure 7r c (black) and the ion pressure tt; (grey). The third column shows the electron 
temperature j? e = Tir a (black) and the ion temperature i9; = Ztiti (grey). The fourth column shows the local values of the normalised 
bremsstrahlung cooling rate (black), the normalised cyclotron cooling rate (grey), and the electron-ion exchange rate (dotted). These 
flows are quasi-one-temperature, and the two-temperature case with small e s are practically indistinguishable from the one-temperature 
flows. 



1/x of the cooling term, then equations l|15|l and 1101 are 
identical for 7 = 5/3, a B = Z and tt c — (1 — t)/x- 



3.1 Quasi-one-temperature flows 

When the electron-ion coupling is strong (i.e. the value of 
tp a i is sufficiently large), we expect the flow to be quasi- 
one-temperature. In Figure we show the examples of flows 



which are effectively one-temperature (with a s — Z and 
tpei = 10) and then how the flows deviate from the one- 
temperature limit when we increase the cyclotron cooling 
rate. Cases shown are those of bremsstrahlung-only cooling 
(e s = 0) , equally efficient bremsstrahlung and cyclotron cool- 
ing at the shock (e s = 1), and dominant cyclotron cooling 
throughout much of the post-shock region (e B = 100). 

For flows with small e s , the electron and ion pres- 
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Figure 2. Stationary structures of two-temperature accretion shock. The system parameters are the same as those in Figure IT1 except 
(rr B ,ip e i) = (Z.0.2), i.e. the electron-ion collisional coupling are weaker. 



sures have the same profiles throughout the post-shock re- 
gion. The temperatures of the electrons and ions are in- 
distinguishable, and the flows are one-temperature. When 
e s is large (i.e. very efficient cyclotron cooling), the elec- 
tron and ion pressure have a different profiles and the elec- 
tron temperature deviates from the ion temperature in a 
small region near the shock. The flow velocity also devi- 
ates from one-temperature model near the shock. However, 
the temperatures of the electrons and ions eventually be- 
come the same further downstream, and the flow is ef- 
fectively one-temperature in the base region of the two- 
temperature cases. Bremsstrahlung X-rays are emitted at 
the high-density regions at the bottom, where the flow is 



practically one-temperature. We would therefore expect the 
X-ray properties to be similar to those of the corresponding 
one-temperature cases. The optical/IR radiation from these 
flows would be somewhat different from those of their one- 
temperature correspondences, as cyclotron radiation (which 
occurs in the optical/IR wavelengths and has a substantial 
optical depth), is emitted mainly from the hotter, less dense 
upper region of the post-shock flow. 

3.2 Two-temperature flows 

In general, where electron-ion exchange is inefficient com- 
pared to radiative cooling (i.e. small ip e i), two-temperature 
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Figure 3. Stationary structures of two-temperature accretion shock. The system parameters are the same as those in Figures [Tl and |21 
except ((TsjVei) = (0.2,0.2). 



effects become significant. Two-temperature flows can also 
occur for moderately large i/> e i, if the initial difference be- 
tween the electron and ion pressures at the shock is sub- 
stantial, i.e. small cr s . In Figures |2] EU and 0] we show three 
examples of the two-temperature flows with various combi- 
nations of the values for the parameters a s and if) e i- 

For (a s ,tp c i) = (Z,0.2), the collisional energy exchange 
with the ions does not keep pace with radiative loss of the 
electrons in the flow. Even though the electrons and ions are 
set to have the equal temperatures at the shock, a strong dis- 
equilibrium (T c < Ti) prevails throughout most of the post- 
shock region. Moreover, while the ions are shock-heated, the 
electron temperature does not rise accordingly due to effi- 



cient radiative cooling and thermal decoupling with the ions. 
The difference between the electron and ion temperatures is 
obvious for all three cases (with different e s = 1) that we con- 
sider, and so is the difference between their electron pressure 
and ion pressure. When cyclotron cooling is weak (e s < 1), 
the velocity profiles of the two-temperature flow and one- 
temperature flows are indistinguishable, but when cyclotron 
cooling is sufficiently strong (say e s = 100) the velocity pro- 
file deviates substantially from that of the one-temperature 
model. 

In all cases, the density is high at the base of the post- 
shock region, and collisional energy exchanges are more ef- 
ficient than cyclotron loss. The collisional exchanges tend 
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Figure 4. Stationary structures of two-temperature accretion shock. The system parameters are the same as those in Figures [Tl and |21 
except (cr s , Vei) = (0-2, 1). 



to bring the electrons and ions towards thermal equilibrium 
here. The bremsstrahlung X-rays are emitted copiously from 
the base of the post-shock flow. In a two-temperature flow, 
the electrons and the ions gradually attain thermal equilib- 
rium near the base because the electron-ion collision rates 
increases with density, and the flow become practically one- 
temperature again. In spite of this, the electron temperature 
gradients near the white-dwarf surface differ between the 
one-temperature and two-temperature flow models. These 
differences are, in some cases, sufficient to affect the line 
and continuum X-ray spectra (see 33] and discussions in the 
later sections). 

We note that the situation can be more complicated in 



the regions near the shock. In the prescription that we con- 
sider, the electron-ion exchange depends on the difference 
between the electron and ion temperatures. Thus, the effi- 
ciency of the energy exchanges between the electrons and 
ions is implicitly determined by the parameter a s . More- 
over, the cyclotron cooling rate, which depends mainly on 
the electron temperature and is most effective in the hot re- 
gion near the accretion shock, is also limited by the efficiency 
of electron- ion exchange and hence a s . The differences in 
properties of cyclotron emission for the one-temperature and 
two-temperature flows should be noticeable. We can clearly 
see these effects in Figures |3 and 2] (cf. Figures Q and The 
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Table 1. The total power of bremsstrahlung radiation L br , an d cyclotron radiation L cy (normalised to p a v^/x s , the kinetic energy of 
the pre-shock flow) for different of (<r s , tpeii e s)' l n an y case with e B = 0, i.e. cyclotron cooling absent, we have L^r = 0-5 an d L cy = 
trivially. 
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0.2 


0.1 


0.1 


0.459 


0.041 


0.2 


0.1 


1.0 


0.346 


0.154 


0.2 


0.1 


10.0 


0.216 


0.284 


0.2 


0.1 


100.0 


0.123 


0.377 


0.2 


0.2 


0.1 


0.451 


0.049 


0.2 


0.2 


1.0 


0.330 


0.170 


0.2 


0.2 


10.0 


0.202 


0.298 


0.2 


0.2 


100.0 


0.114 
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0.443 
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0.2 
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0.1 


0.440 
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0.191 
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0.315 
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0.104 


0.396 
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0.1 
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0.486 


0.014 


0.5 


0.1 


1.0 


0.419 


0.081 


0.5 


0.1 


10.0 


0.288 


0.212 


0.5 


0.1 


100.0 


0.171 


0.329 


0.5 


0.2 


0.1 


0.484 


0.016 


0.5 


0.2 


1.0 


0.408 


0.092 


0.5 


0.2 


10.0 


0.273 


0.227 


0.5 


0.2 


100.0 


0.161 


0.339 
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0.5 
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0.260 
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0.5 
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100.0 
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0.297 
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0.441 
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100.0 
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0.5 
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1.0 
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0.305 
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0.434 


0.066 


1.0 


1.0 


10.0 


0.300 


0.200 


1.0 


1.0 


100.0 


0.178 


0.322 



difference are more prominent when the radiative cyclotron 
cooling rate at the shock, e B , is large. 



3.3 Total radiative loss and X-ray luminosity 

In the hydrodynamic formulation that we consider and 
the boundary conditions that we adopt, all the kinetic en- 
ergy of the accreting gas will be liberated via emitting 
bremsstrahlung X-rays and optical/IR cyclotron radiation. 
We now show that our prescription of the cooling function 
ensures energy conservation and hence self-consistency in 
the hydrodynamic calculations. 

The energy conservation requirement that the total 
power radiated from the post-shock region equals the kinetic 
energy of the pre-shock flow, which is |p a «i, if no energy 
is transported across the white-dwarf surface. We can ob- 
tain the total power of radiation directly by integrating the 
total cooling function over the whole post-shock structure. 
In an explicit representation (with normalised density and 
velocity at the shock), 



At A 

o 7-1 
77- 1 
32( 7 - 1) ■ 



1/4 



dr 



7(1 — r) — r 



7-1 



(20) 



Setting 7 = 5/3 yields the the total power radiated via all 
processes L — 1/2, the result expected for exact energy con- 
servation. 

The power of the bremsstrahlung X-rays is 



1/4 



dr 



7(1 -t) 



7-1 



1 + Acy/Abr 



and the power of cyclotron radiation is simply 

L c y = L — Lbr • 



(21) 



(22) 



Values of Lbr and L cy for different choices of the dimen- 
sionless system parameters are shown in Table Q Values for 
representative choices of white dwarf mass, magnetic field 
and specific accretion rate are given in Table [5] 



4 X-RAY SPECTROSCOPY 

We now calculate the X-ray spectra of two-temperature flow 
models and compare them with the spectra obtained by the 
canonical one-temperature flow models. The question that 
we intend to answer is: how different are the spectral proper- 
ties of two models for a given set of mCV parameters? More 
specifically, will the two-temperature flows produce harder 
X-ray spectra than the one-temperature flows? 

We use the hydrodynamic formulation described in the 
previous sections to generate the density and temperature 
structures of the post-shock emission regions, assuming ap- 
propriate mCV system parameters. In the spectral calcu- 
lation s we adopt the same procedures as in lCropper et al] 
( 1999). We divide the post-shock emission region into a num- 
ber of strata. The strata are each assumed to have con- 
stant density and electron temperature, which take the cor- 
responding mean values in the stratum. The radiative pro- 
cesses in the plasma are collisional-ionisation dominated (we 
ignore photo-ionisation) . Each stratum is assumed to be op- 
tically t hin to the keV X-ray lines and continuum (how- 
ever see IWu et al] 12001). and we may use the MEKAL 
optically thin thermal plasma model llMewe et al] Il985l : 
iKaastra fc Mewd fl99^) in XSPEC to calculate the local X- 
ray spectrum. The total spectrum is the sum of the contri- 
butions from all the strata in the post-shock region (with 
the implicit, simplified assumption t hat line transfer a nd 
scattering effects are unimportant, cf. iKuncic et al.l l2005). 

We consider white-dwarf masses of 0.5, 0.7 and 1.0 Mq, 
and magnetic-field strengths of 10, 30 and 50 MG. The spe- 
cific mass accretion rates m is fixed to be 1.0 erg cm -2 s -1 . 
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Table 2. Parameters of represenative accretion shock models, specified in terms of white dwarf mass, magnetic field strength and the 
election/ion pressure ratio at the shock, cr B . We consider cases with a B = 0.2, cr B = Z (T c = T; at the shock), and a one-temperature 
model (IT) throughout the flow <Wulll994l) . The following columns present corresponding values of the parameters i/> e j and e s , the 
bremsstrahlung and cyclotron luminosities (in units of mVg/x s ), the shock height x s , and the electron temperature at the shock, T o s . 
The plasma composition is approximately solar and the mass flux is set to rh = 1 g cm~ 2 s _1 and the stream has cross-section 10 15 cm 2 . 
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The electron and ion temperatures are set to equal at the 
shock, i.e. a s = Z, and a solar composition is used in deter- 
mining Z and in generating the MEKAL spectra in XSPEC. 

Given the complexity of the lines and the power-law 
like continuum, the spectra of the two-temperature flows 
are not always easy to distinguish visually from those of 
the one-temperature flows, We therefore consider the quo- 
tient spectra, which are the ratios of the spactra of the 
one-temperature flows to the spectra of the two-temperature 
flows. 

Figure |5] shows the quotient spectra (in the 
0.1—10.0 keV band) for various combinations of system 
parameters. These spectra demonstrate that the competi- 
tion between the radiative cooling of the electrons and the 
electron-ion energy exchange affects the post-shock struc- 
ture enough to influence the spectral properties, and the ef- 
fects are quite significant in some cases. Generally, the one- 
temperature model predicts a softer X-ray spectrum for a 
given white dwarf mass. It also produces stronger emission 
lines, especially the Fe L lines. These effects are more signif- 
icant for greater white-dwarf masses. 

The softer X-ray emission of a one-temperature flow 
could be due to the fact that the one-temperature flow tends 
to have a steep velocity gradient in the mid-section of the 



post-shock region, leading to a rapid increase in the elec- 
tron and ion densities. As the X-ray emissivities are propor- 
tional to the square of matter density (for neutral plasmas) , 
this increases the relative contribution of the base region to 
the total the emission. The more realistic two-temperature 
model predicts a more gentle velocity gradient and hence 
density gradient, and substantial X-ray lines and continuum 
are therefore emitted in the hotter strata of the post-shock 
region. When compared with a two-temperature flow, the ef- 
fective temperature of an one-temperature flow is biased to- 
ward the cooler temperatures, and as a consequence a softer 
spectrum results. The excess of the Fe L lines in the one- 
temperature flow can also be explained in the same manner. 

These differences imply that if we fit the observed spec- 
trum of a mCV with both models assuming the same mag- 
netic field and spec ific mass accretion r ate, then the one- 
temperature model JCropper et al.lll999Tl will give a higher 
white-dwarf mass than the two-temperature model. Thus, 
using the one-temperature model will over-estimate the 
white-dwarf mass, and the bias is more severe for systems 
with a strongly magnetised, massive white dwarf. (We defer 
the estimation of observed systems' masses for future work; 
at some level of detail the results may be sensitive to the 
power-law approximation for cyclotron cooling, § 2.) 
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Figure 5. The quotient X-ray spectra in the 0.1 — 10 keV band for various mCV parameters. Each panel presents the ratio of the 
spectrum of the one-temperature flow model to the spectrum of the corresponding two-temperature flow model. In all cases we set 
a s = Z and assume a solar metallicity. The specific mass accretion rate is fixed to be 1.0 g cm" 2 s -1 . The white-dwarf masses are 0.5, 0.7 
and 1.0 Mq (from top row to bottom row). The white-dwarf magnetic fields are 10, 30 and 50 MG (from left column to right column). 
In all cases, the two-temperature flow predicts a harder spectrum than a one-temperature flow. The hardness is more severe when the 
white-dwarf mass and magnetic field increase. 



Now an important question is: how robust are the re- 
sults given that we have assumed a specific <r s ? While the 
value ipei is determined by Z and M w (see Appendix lA20 , we 
treat a s as a free parameter. We carry out calculations for 
two-temperature models with a s — 0.2 and find that they 
produce spectra which differ by less than 2% from those 
where a s — Z (see Figure 

This result can be understood as follows: although a s 
can force the electrons and ions to have unequal tempera- 
tures at the shock, it has weak effects on the flow in the base 
region. This can be seen by comparing the structure pro- 
files in the corresponding panels of Figures [2] and which 
represent models that differ only in their values of a s . The 
density and temperature profiles approach the same values, 
with nearly the same spatial gradients, near the white-dwarf 
surface. Thus, the flow structure at the base is insensitive to 
the electron-ion pressure ratio at the shock. As the majority 
of the X-ray lines and continuum originate from the bottom 
region near the white-dwarf surface, the X-ray properties of 



the two cases should be very similar. Although we cannot 
use the observed X-ray spectrum to infer the value of er s at 
the shock, we can be certain that the mass estimate depends 
on the effective shock temperature and is less affected by the 
choice of cr s in the spectral modelling. 

Inspection of Figure reveals a 'downward' spike near 
the energy of Fe Ka lines, which we identify as the emission 
from the H-like ions. This implies that the two-temperature 
flows and one-temperature flows predict very different ra- 
tios of the H-like (6.97 keV), He-like (6.7 keV) and neutral 
Fe Ka (6.4 keV) lines. The H-like Fe Ka line is weaker for 
the two-temperature flows than for the corresponding one- 
temperature flows, in the systems with a strongly magne- 
tised white dwarf. This effect is stronger for more massive 
white dwarfs. The emission of H-like Fe Ka line requires 
a high plasma temperature (~ 10 keV). The H-like Fe Ka 
lines is expected to originate from regions closer to the shock, 
which have higher temperatures than the region that pro- 
duces the He-like and neutral Fe Ka lines (see e.g. IWu et aTl 
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Figure 6. The quotient X-ray spectra in the 0.1 — 10 keV band: the ratio of the two-temperature models with <r s = 0.2 to that of cr B = Z. 
The mCV parameters are equivalent to those in Figure Only the hard part of the spectrum is affected appreciably, with the case of 
<j b = 0.2 yielding the harder results. The difference at 10 keV is at a level of less than 2% for the case with greatest white dwarf mass 
and magnetic field intensity. 



l200lT) . Thus, the H-like line does not share the characteris- 
tics of the lower ionised Fe lines. Nor, for the same reasons, 
does it share the characteristics of lines of the lighter met- 
als, such as Si, S, Ar and O in the keV spectrum. Moreover, 
as the emission region of the Fe Ka line is relatively near 
to the shock (where the difference between the electron and 
ion temperatures is greatest), the properties of the lines are 
sensitive to the assumed value of a s in the model (see Fig- 
ure 



5 CONCLUSIONS 

We have presented a two-temperature hydrodynamical for- 
mulation for accretion flows in mCVs and used it to cal- 
culate the temperature and density structures of the post- 
shock emission region. Two-temperature effects are signifi- 
cant when the radiative loss is very efficient such that elec- 
trons cannot acquire energy fast enough from the ions via 
collisions. We expect two-temperature flows to occur in sys- 
tems in which the white dwarf is massive and has a strong 



magnetic field. The magnetic field has less influence on the 
two-temperature effects if the white dwarf mass is small, but 
a relatively large effect for more massive white dwarfs. (For 
example, contrast the flow structures in Figures 1 and 2; 
or in Figure 5 consider the greater sensitivity of the spec- 
tra in the bottom row where M = TOM©). In high-mass, 
strong-field systems, cyclotron cooling is more efficient than 
bremsstrahlung cooling. In all cases, the flows eventually 
become one-temperature near the base of the post-shock re- 
gion, where most of the X-rays are emitted. 

In spite of this, two-temperature flows and one- 
temperature flows have distinguishable X-ray properties be- 
cause of the differences in the density and temperature gra- 
dients between the two flows. Our calculations show that 
the X-ray spectra of one-temperature flows are softer than 
the two-temperature flows, if we assume the same system 
parameters, such as white dwarf mass, magnetic field and 
specific mass accretion rate. This result is insensitive to the 
initial difference between the electron and ion temperatures 
at the shock. White-dwarf masses of mCVs obtained by fit- 
ting the X-ray spectra using a one-temperature flow model 
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will lead to overestimates, especially for the systems with 
high white-dwarf masses and magnetic fields. 



ACKNOWLEDGMENTS 

KW and CJS thank Klaus Beuermann for drawing our at- 
tention to the significance of the functional form of the local 
cyclotron cooling. 



REFERENCES 

Aizu K., 1973, Prog. Theor. Phys, 49, 1184 

Beuermann K., 2004, in Vrielmann S., Cropper M., eds, 

ASP Conf. Ser. Vol. 315, Magnetic cataclysmic variables. 

p. 187 

Canalle J. B. G., Saxton C. J., Wu K., Cropper M., Ramsay 

G., 2005, A&A, in press 
Chevalier R. A., Imamura J. N., 1982, ApJ, 261, 543 
Cropper M., 1985, MNRAS, 212, 709 
Cropper M., 1990, Space Sci. Rev., 54, 195 
Cropper M., Ramsay G., Wu K., 1998, MNRAS, 293, 222 
Cropper M., Wu K., Ramsay G., Kocabiyik A., 1999, MN- 
RAS, 306, 684 
Ezuka H., Ishida M., 1999, ApJS, 120, 277 
Fischer A., Beuermann K., 2001, A&A, 373, 211 
Fujimoto R., Ishida M., 1997, ApJ, 474, 774 
Imamura J. N., 1981, PhD thesis, Indiana University, 
U.S.A. 

Imamura J. N., Aboasha A., Wolff M. T., Kent K. S., 1996, 
ApJ, 458, 327 

Imamura J. N., Durisen R. H., Lamb D. Q., Weast G. J., 

1987, ApJ, 313, 298 
Ishida M., Silber A., Bradt H. V., Remillard R. A., Mak- 

ishima K., Ohashi T., 1991, ApJ, 367, 270 
Kaastra J. S., Mewe R., 1993, A&AS, 97, 443 
Kuncic Z., Wu K., Cullen J. G., 2005, PASA, 22, 56 
Lamb D. Q., Masters A. R., 1979, ApJ, 234, L117 
Langer S. H., Chanmugam G., Shaviv G., 1982, ApJ, 258, 

289 

Mewe R., Gronenschild E. H. B. M., van den Oord G. H. J., 

1985, A&AS, 62, 197 
Mukai K., 2003, Adv. Space Rev., 32, 2067 
Nauenberg M., 1972, ApJ, 175, 417 

Potter S. B., Hakala P. J., Cropper M., 1998, MNRAS, 297, 
1261 

Ramsay G., 2000, MNRAS, 314, 403 

Ramsay G., Cropper M., Hellier C, Wu K., 1998, MNRAS, 
297, 1269 

Rothschild R. E., Gruber D. E., Knight F. K., Matteson 
J. L., Nolan P. L., Swank J. H., Holt S. S., Serlemitsos 
P. J., Mason K. O., Tuohy I. R., 1981, ApJ, 250, 723 

Rybicki G. B., Lightman A. P., 1979, Radiative Processes 
in Astrophysics. Wiley, New York 

Saxton C. J., 1999, PhD thesis, University of Sydney, Aus- 
tralia 

Saxton C. J., Wu K., 1999, MNRAS, 310, 677 
Saxton C. J., Wu K., 2001, MNRAS, 324, 659 
Spitzer L., 1962, Physics of Fully Ionized Gases, 2nd Ed.. 
Interscience, New York 



Tennant A. F., Wu K., O'dell S. L., Weisskopf M. C, 1998, 
PASA, 15, 339 

Wada T., Shimizu A., Suzuki M., Kato M., Hoshi R., 1980, 

Progress of Theoretical Physics, 64, 1986 
Warner B., 1995, Cataclysmic Variable Stars. Cambridge 

University Press, Cambridge 
Wickramasinghe D. T., Ferrario L., 1988, ApJ, 334, 412 
Woelk U., Beuermann K., 1996, A&A, 306, 232 
Wu K., 1994, PASA, 11, 61 
Wu K., 2000, Space Sci. Rev., 93, 611 
Wu K., Chanmugam G., Shaviv G., 1994, ApJ, 426, 664 
Wu K., Chanmugam G., Shaviv G., 1995, ApJ, 455, 260 
Wu K., Cropper M., 2001, MNRAS, 326, 686 
Wu K., Cropper M., Ramsay G., 2001, MNRAS, 327, 208 



APPENDIX A: SYSTEM PARAMETERS 
Al Hydrodynamic variables 

The gas is composed of electrons with mass m c plus species 
of ions with masses m\ and charges Z\. If the ionic species 
have fractional abundances f \ by particle number then the 
weighted mean ionic mass, charge and squared charge are 

m = J2i fi m h Z = J2i fi z i and Z* = Ei fi Z i- For example, 
in the specific case of a completely ionised, purely hydrogen 
plasma these constants are m = m p , Z = 1 and Z 2 = 1. 
In these terms, the number density of electrons is 



Z 



m/m c + Z J Tnc 
and the number density of each ionic species is 



Tli= 



fi 



m/m e + Z J m e Z 



9 = tn e 



(Al) 



(A2) 



ensuring a balance of electric charge, n c = ^2 n\Z\. 

For a two-temperature shock, the electron and ion pres- 
sures are unequal. The thermal variables, #{ e ,i}, are given by 



kT a 



m c c 2 n c m c c 2 pc 



P^f m/m^+Z ] (A3) 



kTj 
rmc 2 



Pi 



Pi ( m/m c + Z 
pc 2 \ firm/me 



(A4) 



At the shock surface we define a parameter for the ratio of 
electron and ion partial pressures, a B = P e ,s/ P\,s, with 
summation over ion species i. For a strong shock, the total 
post-shock pressure equals P = P c + Y] jPj = fpavl- As 
the partial pressures are P c = n c kT c and P ; = nifcTi, the 
temperatures of electrons and ions are related at the shock, 



^c,sT c s — (J s ^ ^ TLi s Ti s . 

The post-shock electron temperature is 



3 m e Vf[ 
16 kZ 



/ <?s \ Z + m /m c 
V a s + 1 / T 



(A5) 



(A6) 



and the ion temperature (assumed to be the same for all i) 

(A7) 



T- - —T 
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The upstream pre-shock velocity is assumed to be the 
free-fall velocity at the white-dwa rf surface. We us e the 
white dwarf mass-radius relation of iNauenberd l)l972f) . The 
mass flux of the accretion column m gives the pre-shock 
density p a = rh/va. The shock height x s is calculated by 
equating the bremsstrahlung cooling function at the shock 
with the realistic bremsstrahlung luminosity (involving p a 
and Vff), and substituting the numerically-determined nor- 
malisation constant ip c . 



A2 Electron-ion exchange efficiency ip c i 

When (lAlft . IA2L IA3I and 1A4H are substituted into @ 
and the T e i functions are summed over ion species i then it 
can be shown that the total electron-ion energy exchange is 

r = x P 5/2 p c - 3/2 {p- x Pc) 



(A8) 



where we define X— (Z + 1)/Z and 



X 



W2ive 



lnC* 



m/m c + Z 



5/2 



lnC* 



9.65 x 10 18 — = = Z 2 / 



(m c Z 2 /m) 



5/2 



(A9) 



15 \m/m c + Z y 

in c.g.s. units. Equivalent substitutions in JHJ and summa- 
tion over ion species i leads to the total bremsstrahlung cool- 
ing function 

A br = Ap 2 (P c /p) 1/2 (A10) 
where 



16 



ZP_ ( z 

m e 3 c?h 9B Z \m/m c + Z 



3/2 



5.61 x 10 lo g B Z 2 Z 



—i/2 ( m + Zm c 



-3/2 



(All) 



in c.g.s units. 

In a pure hydrogen plasma with equal elec- 
tron and ion partial pressures (P c = P/2), the 
bremsstrahlung cooling rate is Ab r ~ 3.97 x 
10 16 erg cm~ 3 s _1 (p/g cnT 3 ) 3/2 (P/dyn cm~ 2 ) 1/2 . For solar 
plasma composition (i.e. Z — Zq), we have Z— 1 = 0.09987, 
'Z 2 = 1.391, Z 2 /m = 6.007 x 10 23 g~\ m/m c = 2366. It 
follows that A = 6.99 x 10 16 and X = 2.70 x 10 34 in c.g.s 
units. Thus the general effect of increasing metallicity is 
to increase the efficiency of bremsstrahlung cooling, which 
reduces the shock height if all else is equal. 

The unitless form of the energy exchange function is 

(1 -T- XTTe) 



(7-1) 



(7 - l)V>cV>ei- 



-5/2- 3/2 



(A12) 



and the equivalent unitless function for the bremsstrahlung 
cooling is 



Abr = (7 - 1) V A b>- = (7 - l^cT 



-3/2 1/2 



(A13) 



The dimensionless parameters ip c and ipd (defined in 
llmamura et al.lll99rJ) are constants of each accreting white 
dwarf system in its particular accretion state. Their values 
are, in terms of the characteristics of the accretion flow and 
universal physical constants, given by: 
-4 



(A14) 



and 

lj) c = AlsPa.Vg' 

implying that 

Ai = A^ 2 ' 



(A15) 



(A16) 



which is purely a function of the white dwarf mass and 
composition of the accreting gas. The model parameter ip^ 
roughly describes the rate of electron-ion energy exchange 
compared to the rate of radiative cooling. Upon substitution 
of 1A9II and IA11I we have 



71.2 f\nC\ (Zy m p -= j- 

V 15 J Z 2 m + Zm c . P ' 



9b"8 



(A17) 



where v$ = iiff/10 8 cm s _1 . For accreting magnetic white 
dwarfs lnC w 15 typically, and we take the Gaunt factor 
as g B w 1.25. The factors that depend on gas composition 
tend to lower the value of ip e i for accretion flows of higher 
metallicity. Solar abundances imply a value of ip e i roughly 
32% smaller than for pure hydrogen. 

Using our assumption of solar abundances, the free-fall 
velocity must be greater than vg > 6.21 x 10 s cm s _1 in order 
to make ipa < 1, i.e. radiative cooling comparable to, or 
more efficient than the exchange of thermal energy between 
electrons and ions. For white dwarfs with masses > O.9M , 
the values of ip e i are below 1, and become as low as 0.3 for 
very massive white dwarfs with vs = 1.13 x 10 9 cm s _1 . 
To obtain an extreme value of ip e i =0.1, a free-fall velocity 
of Vfi — 1.96 x 10 9 cm s _1 is required. In this paper we 
consider values of orders tp e i ~ 0.1 — 10, as calculated in 
Figure | A1 | using the white dwarf mass-radius relation of 
3er^|_ 



|Nmi«nb^^|_^^72j) and m ean molecular weight /i w = 2.00. 
In TlmimiuraetHin 1^996), values of ip c i ranging from 0.01 to 
1.0 were, however, used. We note that for white dwarfs in 
mCVs, values of rp e i much greater than 0.01 are necessary. 



A3 Cyclotron cooling efficiency e s 

The timescales of the bremsstrahlung and cyclotron cooling 
processes are respectively 
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It is useful to compare the cooling timescales to express 
the local efficiency of cyclotron cooling with respect to 
bremsstrahlung cooling. The ratio of timescales at any given 
position in the post-shock flow will be written as 

A cv 



(A20) 



This leads to the construction of e s , the relative effi- 
ciency of the cyclotron cooling as evaluated at the shock 
surface, which is an important dimensionless physical pa- 
rameter used for the description and analysis of cooling ac- 
cretion flows in which both bremsstrahlung and cyclotron 
processes are present. Assuming that the cross-section of the 
flow is circular, e a is obtained in realistic terms by substi- 
tuting appropriate system parameters and shock conditions 
into the equations for the A cooling functions and 1A20II . 



14 C. J. Saxton et al. 



0.2 0.4 



0.6 0.8 1.0 

M, / Me 



1.2 1.4 



10* - 



10 2 - 



10° - 



1 0" z - 5MG 




10- 



0.2 0.4 



0.6 



0.8 
/ Mo 



1.0 1.2 



Figure Al. Electron-ion exchange efficiency parameter ip e i 
verses white-dwarf mass. We assume approximately solar abun- 
dances. This parameter is independent of the magnetic field 
strength and a s . 



In general, the efficiency of cyclotron cooling relative to 
bremsstrahlung cooling, in terms of electron temperatures 
and number densities, is 



Figure A2. Cyclotron efficiency parameter e s as a function of 
white-dwarf mass, given for five choices of magnetic field strength. 
Here we set the ratio of pressures at the shock as cr s = Z; the 
mass accretion rate is rh = 1 g cm~ 2 s^ 1 ; the pole area is a = 
1.0 X 10 15 cm 2 . The flow cross-section is assumed to be circular, 
and the composition is solar. 



e = 0.0762- 



,Z 2 



5 D 2-85 
i>7 n cl6 



(A21) 



Since a circle is the two-dimensional geometric shape with 
the minimum ratio of perimeter to internal area, the above 
expression for e is actually a lower limit. For more realis- 
tic cross-sections, the numerical factor 0.0762 would become 
larger. 

By substituting expressions relating n e i6 and tempera- 
ture 1A6II to the pre-shock density p a and free-fall velocity 
vg, the efficiency e a can be re-expressed in terms of the prop- 
erties of the white-dwarf accretion flow. In the general case 
with multiple species of ions 

2.13 x 10~ 1(i (Z + m/m B ) 3 - 85 
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where p a -8 = p a /10~ 8 g ■ cm~ 3 and vg = 
This expression is slightly different 
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v tt /10 8 cm ■ 
from the 



cy- 

clotron /brem s strahl ung efficiency ratio given in 
lLanger et alJ il982f) where the emission region is semi- 
infinite but without a specific geometry. Assuming a 
circular flow cross-section and parameters appropriate for 
accretion shocks in AM Herculis systems, with 015 = 1 
and m = 1 g cm _2 s _1 , the cyclotron-cooling efficiency 
parameter varies with white-dwarf mass and magnetic field 
as shown in Figures IA2I and IA3I 
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Figure A3. Same as Figure IX2I but with <r s = 0.2. 



Two-temperature accretion flows in mCVs 15 



A4 System geometry 

The geometry of the accretion stream could be important 
as optically thick cyclotron cooling depends on the emitting 
surface area. Under our approximation, 1A21L the efficiency 
parameter e s is proportional to the ratio of perimeter to 
cross-sectional area. 

A circular cross section for the accretion column gives 
the minimum value of e s - All else being equal, flows 
with cross-sections departing from a circle cannot have 
lower values of e s than given by 1A22II . The greater the 
non-circularity, the greater e s must be. There is observa- 
tional evidence that in AM Herculis type systems (po- 
lars) the cross-section of the accretion col umn near the 
the white-dwarf surface is a banana-like arc dCropperlll985l 
IWickramasinghe fc Ferraridlf 98sl ; iPotter et al]ll998l) . This 
geometry ought to give a e s value several times higher. 
Whether the real banana-shaped cross-sections have more 
crenulated or convoluted edges at finer spatial scales is un- 
known. If they do, then the equivalent e s values may be 
considerably greater than the lower limit of the circular ap- 
proximation. Because of the effect of the surface area to 
volume ratio in a system radiating via an optically thick 
process, real effective values of e s may be less sensitive to 
the magnitude of the polar magnetic field than to other in- 
fluences determining the flow structure, e.g. the geometric 
relationship between the magnetosphere and the companion 
star, and the condition of the flow where it threads onto 
the magnetic field. If flow geometries are sufficiently diverse 
among mCV systems then there may exist some strong-field 
systems which have lower e s than weaker-field systems that 
have more circular flow sections. The magnetic field strength 
is measurable, implying a lower limit on e s in each system, 
but unfortunately the value of e s is not a directly observable 
quantity. 



